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ABSTRACT 

We examine the dust distribution around a sample of 70,000 low redshift galaxy groups and 
clusters derived from the Sloan Digital Sky Survey. By correlating spectroscopically identi- 
fied background quasars with the galaxy groups we obtain the relative colour excess due to 
dust reddening. We present a significant detection of dust out to a clustercentric distance of 30 
Mpc/h in all four independent SDSS colours, consistent with the expectations of weak lensing 
masses of similar mass halos and excess galaxy counts. The wavelength dependence of this 
colour excess is consistent with the expectations of a Milky Way dust law with Ry=3.1. Fur- 
ther, we find that the halo mass dependence of the dust content is much smaller than would be 
expected by a simple scaling, implying that the dust-to-gas ratio of the most massive clusters 
(~ 10" h-^ Mq) is - 3% of the local ISM value, while in small groups (- lO^^-^ /i^^ Mq) 
it is ^ 55% of the local ISM value. We also find that the dust must have a covering fraction 
on the order of 10 % to explain the observed color differences, which means the dust is not 
just confined to the most massive galaxies. Comparing the dust profile with the excess galaxy 
profile, we find that the implied dust-to-galaxy ratio falls significantly towards the group or 
cluster center. This has a significant halo mass dependence, such that the more massive groups 
and clusters show a stronger reduction. This suggests that either dust is destroyed by thermal 
sputtering of the dust grains by the hot, dense gas or the intrinsic dust production is reduced 
in these galaxies. 
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Dust grains have long been known to play an important role in the 
interstellar medium and the star formation which occurs within a 
galaxy. Further, because dust grains can absorb and redden back- 
ground sources, an accurate knowledge of the large scale distri- 
bution of these grains is crucial. Groups and clusters of galaxies 
present a unique opportunity to study this large scale distribution as 
well as the processes important in dust evolution. Dust within mas- 
sive clusters is thought to sputter on timescales of 10^ - 10^ years 
jDraine & Salpe ter 1979), depending on the density and tempera- 
tures of the environment. This is caused by the ejection of atoms 
from the dust grain by the collision with sufficiently energetic gas 
particles. Given the presence of a distributed, hot plasma in groups 
and clusters, the short timescale for sputtering means that dust ob- 
served in these systems must have been accreted recently. This po- 
tentially gives a probe of the relevan t dust creation processes. 

Recentlv lMenard et al.l ( l2009bl) have shown that dust excesses 
exist out to Mpc scales around i<21 galaxies. There are a vari- 
ety of mechanisms by which dust may escape from galaxies and 
become distributed on such large scales. Star-forming galaxies are 
often seen to have large outflows of gas and dust which are blown 
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out by the power fro m supernova feedback dHeckmanetalll 19901 : 
iTremonti et alj l2007h . In addition, active galactic nuclei (AGN) 
can have significant power, and the jets they induce may be able 
to remove gas and dust fro m galaxies and redistribute it within 
the larger environment (e.g. iMcNamara & Nulsenll2007h . During 
galaxy-galaxy mergers a significant amount of the gas and dust may 
be removed due to collisional processes, or su bsequently blown ou t 
by the induced star-formation or AGN power ('Hopkins et al.'2006). 
Observations of galaxies falling into massive clusters show gas and 
dust being stripped fro m the galaxy and incorporated into the su r- 
rounding environments jCrowl et al.ll2005l : lDomainko et al.ll2006h . 

The temperatures and densities of the large scale environments 
of groups and clusters are difficult to probe observationally. While 
X-ray observatories have allowed the determination of the tem- 
peratures and densities o f massive clusters to near the virial ra- 
dius (e.g. ivikhlinin et allboosh . analysis of the detailed properties 
of representati ve samples of gr oups has largely relied on stacking 



large samples jDai et al" 



20071) and/or bee n limited to measuring 



total X-ray luminosities i Rvkoff et alj2008l) . As a result, the densi- 
ties and temperatures of the environments o utside the virial radius 
are largely infer red from simulations alone jPfrommer e7Zll2006l : 
iKav et al.[|2007l) . Indeed, the non-detection of significant amounts 
of baryonic material in the local universe has given rise to the pos- 
tulate that this material is contained in overdense gas with tem- 
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peratures of 10^-10^, dubbed the warm-hot intergalactic medium 
(WHIM) (see lBregmarJl2007h . The expected temperature and den- 
sity of the WHIM are likely able to destroy dust through thermal 
sputtering. Therefore, measuring the large scale dust distribution 
is a potentially powerful probe of the large scale temperatures and 
densities. 

Previous attempts to measure the dust centered on clusters 
has focused on the inner regions, largely within the virial radius. 
Observations of typical groups (< 10^"^'^ M©) have not yet 
been done. The first attempts to quantify the dust content of clus- 
ters involved the counting of relative deficiency of background 
sources. These attempts general ly agreed that the 1/— band atten- 
uation is on the order of 0.2-0.4 fewickv" 1 95 1 : Bo gart & Wagoned 
Il973l : iBovle et alj|l988l ; iRomani & Maoz 1992) . However, early 
infrared observations of clusters attempting to detect direct emis- 
sion from this dust largely led to non or marginal detections, im- 
plying dust masses much lo wer than that implied by background 
counts dAnnis & Jewitt 19931. IWise et al.ll 19931 : [Stickel et ali2002l : 
iMontier & Giard.2005) . Similarly, recent observations with Spitzer 
have also faile d to show dust m asses implied by 0.2-0.4 magnitudes 
of attenuation jSai et al.ll2007l) . 

Recently, large, uniform surveys have made detecting statisti- 
cal colour ex cesses of backgro u nd sou rces of large samples of clus- 
ters feasible. IChelouche et alj j2007h have correlated background 
quasars with intervening galaxy clusters and found that, for mas- 
sive clusters, the excess reddening, E(B-V), is on th e order of 
0.004 magnitudes in the central Mpc. IBovv et al.l fcoOSh . using the 
spectra of early type galaxies, found an upper limit on the extinc- 
tion within 2 Mpc of massive clusters t o be Ay < 0.003 m agni- 
tudes. Similar results were obtained by iMuUer et alj ilOOA who 
found, using photometric redshifts to identify background galaxies, 
a similar result of Ay = 0.004 ± 0.010 mag of attenuation. How- 
ever, these type of measurements are differential measurements 
and, as such, are dependent on the control sample of sources with 
which the colour excess is measured. Since it is known that clus- 
ters of this ha lo mass have s ignificant excess ma ss out to at least 
20 Mpc dSheldon et^[2007b ). and indeed IChelouche et al.l 
( I2OO7I) has shown that the colour excess varies at least out to 5 Mpc 
from the cluster center, then the presence of dust on these larger 
scales may bias the measurement of the central regions. 

In this paper, by measuring the reddening effect on back- 
ground quasars, we examine the large scale (~ 50 h^^ Mpc ) 
radial profile of dust centered on groups and clusters. This is the 
first measurement of such large scales and the first to probe such 
low mass groups. The large scale is important to separate the clus- 
ter/group dust from the dust expected to be associated with indi- 
vidual galaxies. We discuss the data in ij2]and present the method 
and the measurement in ij3] We discuss the implications of the mea- 
surement in S|4]and summarize our findings in S|5] Throughout this 
paper, we adopt, as was done during the assembly of the group cat- 
alogue, a ACDM cosmology with the parameters of the third year 
WMAP data, namely = 0.238, = 0.762, fib = 0.042, 
n = 0.951, h = ffo /(100kms~^ Mpc~^) =0.73 and ag = 0.75 
dSpergel et alj2007l) . 



2 DATA 

The technique we use to examine the dust content of galaxy groups 
and clusters is conceptually simple. Our approach relies on mea- 
suring the change in the mean colour of background sources as a 
function of their projected distance from groups and clusters. It is 



expected that the change in the mean colour will be small, and thus 
we must stack the signal from many clusters together to measure 
the effect on the distribution of background sources. Therefore, we 
require a large sample of groups and clusters, which are uniformly 
selected and have well defined masses. We also require that we have 
a large number of background sources with an intrinsic colour dis- 
tribution with little scatter and which do not vary with sky position. 
The background sources must also be at high redshift, so that they 
are not physically associated with the clusters we are examining. 
For these purposes, the best publicly available data is derived from 
the the Fourth Data Release (DR4) of the Sloan Digital Sky Survey, 
a five colour (ugriz) photometric and spectroscopic survey cover- 
ing over 4780 deg ^ and containin g ^ 670,000 spectra of galaxies, 
quasars and stars jAdelman-McCar thv et al. 2006). Below we dis- 
cuss the group sample and the sample of background objects drawn 
from this survey which are used in this paper. 



2.1 Galaxy group catalogue 

The majority of galaxies in the local univers e reside in some kind of 
association with at least one other galaxy ( iEke et af]|2004h . How- 
ever, quantifying the mass of a large sample of those associations 
is difficult. In particular. X-ray temperatur e or luminosity, often 
used as a mass indicator in galaxy clusters dReiprich & B6hringe3 
l2002h . is too low or faint in galaxy groups to determine mass for 
a survey the size of the SDSS. Also, the velocity dispersions of 
galaxies, which are often used as tracers of the potential well of 
clusters or massive groups, is ineffective when the groups contain 
only a handful of spectroscopically confirmed galaxies. Recently, 
lYangetaljj200^ showed that an effective mass estimate of galaxy 
groups can be obtained by essentially ranking groups by their total 
luminosity and associating these rankings with the expectations of 
a A CDM halo o c cupat ion model. 

lYang et alj ( |2007|) have applied their algorithm to the SDSS 
DR4 to produce a sample of ~ 300,000 galaxy groups with masses 
as low as 10^^'^ Mq. In this paper we use "Sample I" from 
lYang et alj ilOOf). which exclusively uses galaxies with SDSS 
spectroscopic redshifts. We chose this sample, as opposed to the 
other samples, which add in existing redshifts from the literature, 
because we are principally concerned with obtaining a uniformly 
selected population of galaxy groups. In the most recent galaxy 
group catalogue of lYang et a l.. the authors rank the galaxy groups 
both by total luminosity in the r— band as well as the total stel- 
lar mass of the group. In this paper, we will use the group masses 
obtained by the ranking the total stellar mass, with the goal of min- 
imizing the effect of a particular group's recent star formation his- 
tory. However, we note that because we must stack many groups 
in relatively wide mass bins, this choice has no effect on the re- 
sults. Finally, we restrict our analysis to galaxy groups with masses 
greater than 10^^'^ M0, which leaves a final sample of 75947 
groups between z=0 and z=0.2. Figure [T] shows the halo mass dis- 
tribution of our group sample broken into three mass bins which 
we use later in the paper. 

2.2 Background objects: spectroscopically identified quasars 

Our principal concerns when choosing a sample of background 
objects to examine are that the sources are at sufficiently high 
redshift that they are not physically associated with the group, 
that they represent a relatively uniform population, with a small 
dispersion about their average properties, and that their colours 



Dust accretion and destruction in galaxy groups and clusters 3 




1013 1014 1015 

Halo Mass (Mgh"') 



Figure 1. The halo mass distribution of the galaxy group and cluster sam- 
ple. The sample is shown in three groups of halo mass, from lO^'^'^ to 
10^3 ft-l M0(red), lO" to lO^^ ** h'^ MgCblue), and IQi^-S to IQl^ S 
Mq (black). The short dashed (long dashed) line in each bin represents 
the median (mean) halo mass of the bin. 

are well calibrated. For all these reasons we have chosen to ex- 
amine the photometric properties of spectroscopically identified 
quasars drawn from the fourth edition of the SDSS quasar cata- 
log ( Schneider et al. 2007). This catalog contains ~ 77,000 quasars 
brighter than Mi = —22 drawn from the fifth SDSS data release, 
and are targeted based on their photom etric properties and/or th e 
presence of an unresolved radio source jRichards & others]|2002h . 
We reduce this sample to contain only quasars within DR4, ie. those 
for which we have an identified group sample. The SDSS quasar 
catalog contains pi/ magnitudes for each of the five bands of SDSS 
photometry (ugriz), with typical errors of 0.03 mag. The colours 
have been corrected for galactic extinction using the dust maps of 
ISchlegel et alj ( Il998l) . We only examine quasars with redshifts be- 
tween z = 1 and z = 3 in order to avoid quasars physically associated 
with the target groups, or extreme objects at high redshift. 



3 ANALYSIS 

Our first goal is to measure the change in the mean colour of the 
background quasars as a function of their position from foreground 
clusters. This change in the mean colour is also known as the excess 
reddening because interstellar dust is observed to extinguish light 
more readily at the blue end of visible light, causing the light to 
redden. 

We parametrize the excess reddening, E{i — j), of a source, 

k, as 

E{i - j)'' = rrii - m] - {rm - ■mj)controi (1) 

where and rrij are the magnitudes in the i and j bands of the 
source for which the excess reddening is measured. The final term 
in this expression is the mean colour of the control sample of ob- 
jects. In practise, the excess reddening measured in this way of a 
single background source is dominated by the intrinsic distribution 



width of colours of that source. However, with the assumption that 
the background sources have a well defined mean colour we can 
stack large numbers of sources to measure a mean excess colour, 
{E{i — j)). Therefore, by stacking many backgrounds sources in 
bins of projected distance, r, from the centre of our group and clus- 
ter samples, we can obtain an excess colour profile of the groups 
and clusters, {E[i — j)){r). 

A key feature of this type of measurement is the control sam- 
ple of sources. Unless the unobscured mean colour of the back- 
ground sources is known, it must be made as a differential mea- 
surement with respect to a control sample of backgrou nd sources. 
In the literature, the control sample is not well defined. IBovv et al] 
1 I2OO8I) use a control sample of all background sou rces greater than 
2 Mpc from the center of a cluster, while IChelouche et al.l 
( l2007h defines the control sample to be background sources which 
are greater than IR200 from the cluster center. We have, by trial and 
error, found that the differential dust signal approaches zero only at 
>40 h^^ Mpc from the center of the cluster. So, the control sam- 
ple of background sources is must be taken from large radius. The 
control sample of quasars are defined as those which have projected 
clustercentric distances, r, of 46 h^^ Mpc < r < 50 h^^ Mpc from 
a given cluster in the sample. 

Notice that, similar to the weak lensing or cluster correlation 
measures, our dust measurement is a correlation measurement. A 
given quasar can have a small projected distance from cluster A, 
but still be counted as a control quasar if it is in the control range 
of cluster B. In other words, each quasar is in many clustercentric 
radial bins. Because of the much larger area covered by the con- 
trol sample of quasars, the number of quasars in the control sample 
greatly outnumber the number of quasars in any of the other sam- 
ples. This enables us to do a Monte Carlo estimation of the error in 
the excess colour by selecting 200 independent samples of control 
quasars which have the same number as a given target bin. The er- 
ror bars are then the 1 sigma limits of the excess color from the 200 
trials for each bin. 

3.1 Large scale distribution 

We begin by examining the total colour excess from a stacked set 
of all the groups in our sample. We measure the mean colour of the 
background quasars within clustercentric annuli in projected radius 
bins of 2 h^^ Mpc , out to 50 h^^ Mpc from the cluster center. 
The five band photometry of the SDSS quasar sample allows for the 
measurement of four independent excess reddening signals. Figure 
|2] shows the colour excess for each of the four sets of independent 
colours. 

As Figure |2] shows, a significant colour excess is measured 
in all four independent colours as far out as 30 h^^ Mpc from 
the group and cluster centre. The typical projected virial radius 
of a group or cluster of this mass (M ~ 10^^'^ M©) is only 
~ 1 h^^ Mpc . Thus, such a significant large scale distribution 
of reddening is initially surprising. However, studies of the cluster 
correlation function (Bahcall & Soneira 1 9831), and the weak lens- 
ing profile of similar clusters dSheldon et al.ll20()7bl ; |johnston et al.l 
I2007D shows that there is a significant excess of matter out to sim- 
ilar distances from clusters. We expect that the dust excess to such 
large scales is just a result of this excess matter, although we will 
address this further in i|4.2| 

While the size of the measured excess reddening shown in Fig- 
ure[2]is quite small (< 0.015 mag), if it is due to the presence of dust 
then it corresponds to a large and extended dust distribution. How- 
ever, there are some systematic effects which must be accounted 
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Figure 2. The excess colour of four independent colours as a function of distance from the group center The excess colour is measured with respect to the 
colour of quasars at a projected distance of 46 to 50 Mpc . The error bars are 1 sigma errors from Monte Carlo estimation. The blue, dashed line (red, 
dotted line) represents the Ry=2. 0(5.0) dust reddening law which minimizes the chi-square for Ay as a function of radius. 
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for before we can be sure the reddening is due to dust. In particu- 
lar, galaxy groups and clusters have significant mass, and therefore 
sources behind this mass will be gravitationally magnified. A sam- 
ple of sources which is chosen by a fixed flux limit will lead to 
more high redshift sources, and a lower absolute luminosity limit 
at fixed redshift behind the cluster, than for a patch of sky far from 
the cluster. Thus this measurement is potentially affected when the 
intrinsic mean colour of the background sources is luminosity or 
redshift dependent. 

We can estimate the size of the magnification effect of our 
groups and clusters. We would expect this effect to trace the mass 
distribution, so it wo uld be most pronounced close to the cluster. 
Ijohnston et alj ( l2007h measured the mass profile of a large sample 
of galaxy clusters and found that at a distance of 1 Mpc from 
the most massive clusters the surface mass density is ~ 10^ 
h Mq/pc^. This corresponds to a magnification of ^ 0.02 mag- 
nitudes — a minute change in the effective limiting magnitude of 
our quasar sample. Given this, we see that for our <g-r> mea- 
surement of ~ 0.011 at 1 Mpc to be completely explained by 
magnification of a quasar sample with varying mean colour, the in- 
trinsic cumulative g — r colour of our quasar sample would have 
to vary by 0.55 per magnitude near the magnitude limit. However, 
we find that the cumulative g — r colour of the control sample of 
quasars only varies by 0.005 magnitudes from i=19 to i=20.2, the 
magnitude limit of our sample. In other words, as it is for all four in- 
dependent colours, the size of the magnification-induced reddening 
effect is two orders of magnitude smaller than required to explain 
our results. 

We note that the possibility of having foreground emission 
from the group or cluster contribute to the the reddening signal is 
not physical, given that the local background subtraction used in 
the quasar photometry will remove this component. 



A value of Rv = 3.1 is generally used based on studies of inter- 
stellar Milky Way dust, but values in the r ange from 2.5 to 5 have 
been measured for t he Magellanic cloud s jPrevot etal.]|l984h and 
starbursting galaxies jCalzetti et alj2000l) . In contrast, the emission 
from a typical collection of cluster galaxies would have a negative 
Rv'. Unfortunately, we can not directly determine Rv because we 
do not have B and mag nitudes. However, using the fitting func- 
tions of lO'Donnelll ( Il994l) . we can transform our excess colours to 
determine Ay and to determine whether the shape of the expected 
reddening curve is consistent with the known properties of inter- 
stellar dust. 

In Figure (2] we also show the predictions of two dust redden- 
ing laws, one with Kv =2.0 and one with Rv=5.0. To do this we 
must find the Ay value at each radial bin from the combination 
of the four independent colours. We do this by minimizing the chi 
square value at each step given the errors on each colour For the 
majority of the colours and radii, these two dust laws essentially 
bracket the data. This seems to suggest that the colour profile is 
very similar to that expected from known dust laws. However, we 
would like to make a more precise measurement, which we do by 
stacking the data in radial bins. 

In Figure [3] we show the reddening of background quasars in 
two radial bins for each photometric band with respect to the r band. 
The solid black line in each panel is the extinction curve expected 
from an interstellar dust law with Rv=3.1. The dust attenuation, 
Av, is 0.0223 in the panel showing the inner 10 Mpc , while it 
is 0.0097 for the 10 > r > 20 h'^ Mpc bin. While the Rv=3.1 law 
shows reasonable agreement with the wavelength dependence of 
the reddening, we also allow Ry to be fit simultaneously with Ay . 
These results are shown with the dotted line and correspond to best 
fitting Kv values of 3.3 and 3.5 in the two radial bins respectively. 
Thus, the reddening is very similar to the expectations of a Milky 
Way dust law. 



3.2 Wavelength dependence of reddening 

Here we consider the wavelength dependence of the reddening, to 
gain some indication of the nature of the dust. The wavelength de- 
pendence, Rv is commonly parametrized by linking the absolute 
extinction in the V band, Av to the excess B — V colour as 



3.3 Halo mass dependence 

We now explore how this dust distribution depends on the mass of 
the host halo. In Figure|4]we show the qso <g-i> colour excess as 
a function of distance from the group center in three bins of total 
group mass. The <g-i> colour is chosen as our principle measure 
of the dust reddening in the remainder of the paper. These two pho- 
tometric bands are the best calibrated, and have been recently used 
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Figure 3. The colour excess in two bins of radial distance from the group 
centres. The solid bla ck line is the expe ctations from a Rv=3.1 dust law, us- 
ing the expansion of iO'Donneiilll994l) . with a attenuation of Ay = 0.0223 
and 0.0097 for the > r > 10 h'^ Mpc and 10 > r > 20 h'^ Mpc bins 
respectively. The dotted line is the result of simultaneously fitting for Rv 
and Av, which results in (RvAv) = (3.3, 0.0220) and (3.5, 0.0094) in the 
two radial bins, respectively. The data points are for the colours u-r, g-r, r-r, 
i-r, z-r respectively. 



in the literature. Conveniently, the <g-i> is mo st closely related to 
the more c ommonly used B — V dust colour as jPrevot et al.l 19841 : 
iMenard et al.,.200&) 



E(a -i) = 



1.55E{B ~ V). 



^E{B - V) 



(3) 
(4) 



Using this transformation, and assuming Rv=3A, we show 
the corresponding dust attenuation in the V band, Av, on the right 
ordinate axis. Quite strikingly we see that the centers of the clusters 
have dust attenuations of Av ^ 0.03-0.04, when measured with re- 
spect to quasars ^ 45 Mpc away. In comparison, the disk of a 
spiral gal axy at a similar redshift has a dust attenuation of Av ~ 
0.25 mag jHolwerda et al . 2009). This suggests that the dust signal 
cannot be localized to disks within the line of sight alone, because 
the covering fraction of such group and cluster members is much 
less than (0.03/0.25) = 0.12. In other words, there must be a sig- 
nificant component of dust which is not localized within the disks 
of ~ 0. I L* galaxies alone. A similar conclusion was reached re- 
cently bv lMenard et"ai] j2009bh . who showed that the dust attenua- 
tion around galaxies extends well beyond the radius expected if the 
dust was confined to the disk. We examine this further in H3A\ 

It is also noticeable that the dust distribution shows a relatively 
small halo-mass dependence, in that the most massive bin shows a 
larger colour excess than the lowest mass bin at each position out 
to 20 Mpc . However, except at the inner most bins, the mea- 
sured colour excess in each of the three mass bins is always within 
~ 0.001-0.002 magnitudes, a much smaller range than the range 
of their values as a function of group centric distance. A possible 
explanation for such a small dependence on the halo mass could be 
that the halo masses themselves are uncertain, so that perhaps the 
wide halo mass bins actually contain roughly the same size groups. 
However, it is worth noting that the fraction of red galaxies within 
these same groups show a significant halo mass dependence. 50% 
of the Lt galaxies in 10^* h^^ M0are "early type" while only 
20% of the same galaxi es are "early type" in 10^^'^ MQhaloes 
dWeinmann et al.ll2006l) . This argues that the masses are well de- 
fined for a statistical study such as this. 

In the bottom panel of Figure [4] we show the same <g-i> 
colour excess but now as a function of scaled radius, namely Ri8o- 



|o.os 

E0.015 

u 

v 

I, 0.01 
"3 

"0.005 

A 

M 



"T" 



10'3-8 Mgh-' 
lO'^o Mgh-' 
10'=-^ Mgh-' 



i i 



1 1 



! i 



i i 



! i i i i f I 



i ! •- 



0.04 

0.03 ' 

I 

0.02 , 
0.01 ' 




10 20 30 40 

Distance from group centre (h"' Mpc) 



1 0-02 f- 
S 0.015 



" 0.005 



— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 



» i . 



' ■ 1 • , I I • 

_l I I I I I I I I 



- 0.04 
0.03 
0.02 
0.01 




10 20 30 

Distance from group centre (R,go) 

Figure 4. The g — i colour excess in three samples of galaxy groups as 
a function of distance from the group center. The V band attenuation Av 
is plotted on the right ordinate axis assuming Rv=3.1. This is shown as a 
function of physical distance (top panel) and in Rigo (bottom panel). 



For this plot we have used the individual Ri8o of each cluster, rather 
than using the median Riso of each bin. This is an important point, 
as simply scaling physical distance in the top panel by the median 
Rigo would result in lower mass groups actually have a higher red- 
dening than higher mass groups. However, because there is a range 
of Ri8o in each bin, the result is to almost completely remove the 
halo mass dependence, especially at distances far from the cluster. 
It appears that the dust reddening is similar at a given scaled radius 
far from the cluster. We will return to the interesting behavior at 
small scales in 34.21 



3.4 Spatial distribution of the dust 

By measuring the quasar colour excess as the difference in the mean 
value of the quasars, our results are insensitive to the spatial distri- 
bution of the dust. If the dust exists in concentrated clumps with 
a small covering fraction, then the colour excess could be driven 
by relatively few objects with large reddening values. However, if 
the dust is uniformly distributed we would expect each quasar to be 
reddened by the mean value. 

In an attempt to address this question, we explore the color 
excess as a function of the percentiles of the quasar color distribu- 
tions. In order to determine what fraction of quasars contribute to 
the signal, we sort the quasars in a particular radial bin and those 
in the "control" background sample by their measured g-i colours. 
We bin each sample by percentile and measure the mean g-i colour 
in that percentile bin, and then subtract the measured mean g-i per- 
centile colour of the cluster quasars from the measure mean g-i 
colour of the "control" sample for in the corresponding percentile. 
This leaves us with the {g — i) colour excess in percentile bins. 

Admittedly, this measure is difficult to interpret physically, so 
we first attempt to develop a framework in which to interpret the 
results by using simulated reddening on background quasars. We 
take as our control sample of quasars those which are the control 
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Figure 5. Tlie simulated g — i colour excess as a function of percentile 
for four different reddening models. The model is set to have a mean color 
excess of 0.01 and F^gjj determines the fraction of objects which are red- 
dened. 



sample of the lO^'^ '^-lO^'^ mass bin. We then add an amount of 
reddening to some fraction of the background quasars (Fred) such 
that the mean colour excess is 0.01; for instance, all quasars can be 
reddened by 0.01 magnitudes, or 1/3 of the quasars can be reddened 
by 0.03 magnitudes. The results are shown in Figure[5] where Fred 
is varied to be 1, 0.15, 0.05 and 0.01. As expected, the color excess 
as a function of color percentile is flat for Fred=l and becomes 
more dominated by a high percentile peak with lowered Fred- 

We now present Figure (6] which shows the color excess as 
a function of the percentiles of the real quasar color distributions. 
This is shown for several radial bins and for each halo mass range. 
We notice that for the majority of the radial bins the quasar color 
excess is relatively constant, never more than ~ 50 % away from 
the mean value in that bin. Comparing these with the results of 
our simulated data suggest that 15 % of the quasars are causing 
this signal. Further, we see that the color excess curves become 
progressively steeper with decreasing radius, implying that the bulk 
of the signal is coming from fewer objects as we get near to the 
cluster In the inner-most bin the reddening can be explained with 
5-15% of the quasars dominating the signal. 

As we will see in section ii4.1l there are ~ 20 galaxies (i < 21) 
in the central 2 Mpc of our most massive clusters. Thus, as- 
suming that the reddening signal is concentrated uniformly in "bub- 
bles" around these galaxies, a 10% covering fraction requires the 
"bubbles" to have radii of ~ 140 h'^ kpc . These bubbles are sig- 
nificantly larger than the typical size o f the stellar mass in a galaxy 
dShen et al.ll2003l : lMcGee et al.ll2008l) . So, while this excess dust 
is not isotropically distributed throughout the groups, it also is not 
confined to massive galaxies. 

The results of Figure [6] are also relevant in understanding a 
possible selection effect on our results. Our colour excess measure- 
ment are based on samples from a magnitude limited survey, and 
thus dust extinction may move quasars out of our sample, thereby 
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Figure 7. The excess dust surface mass density, in units of h Mq/pc^, as 
a function of group-centric distance. This is shown for three different bins 
of group mass. 



leaving a relatively reduced color excess. However, examining Fig- 
ure|6]we see that there does not exist any significant population of 
highly extincted quasars. Further, this result exists even when the 
sample is reduced to different magnitude bins. The lack of a highly 
extincted population suggests that this is not a significant effect for 
the bulk properties of the sample. 



4 DISCUSSION 

We have shown that there is a significant distribution of dust cen- 
tered on galaxy groups and clusters and, in this section, we would 
like to measure some of its physical properties. This requires more 
information about the physical nature of the dust and, while we 
were able to show that the reddening is consistent with that ex- 
pected from local interstellar dust, our measurements are somewhat 
crude. In what follows we will therefore make the assumption that 
the dust is similar to Milky Way interstellar dust with an Rv=3.l. 

We are most interested in converting the Ay values we have 
measured into a dust mass surface density. To do this we need an 
estimate of Kext, the extinction cross section for a given mass of 
dust: 

Av 



TjMdust ~ 



Kext 



(5) 



For this , we use the carbonaceous-sili cate dust grain model d evel- 
oped bv lWeingartner & Dra ine' ('200 1*) and lLi & Drain3 j200lh and 
subsequently tweaked by [Draine (2003!)- These models have been 
shown to successfully reproduced the observed infrared emission, 
scattering properties and extinction properties of local interstellar 
dust. The value of Kext is wavelength dependent, but in the V band 



it is iCe^t = 1.54 X IQ'^cm-' 



. With this assumption we can now 



directly show the surface mass density of dust which is responsible 
for the reddening signal. We plot this in Figure [7] for each of the 
three halo mass bins. 

From this we conclude that the dust mass located within 1 
h" Mpc of the cluster center is on the order of lO"' h"^ Mq, 
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Figure 6. The relative color excess at a given percentile of the quasar color distribution. This is shown in several radial bins for each halo mass range. The 
color excess is measured with respect to the color at a given percentile of the control quasars. The red, solid line is for quasars within the central 2 Mpc , 
while the blue, dotted line is for quasars within 9 < r ^ Mpc ) < 1 1 . The green, dashed (black, long dashed) line is for quasars within 19 < r (h~ ^ Mpc ) 
< 21 (29 <r(/i-i Mpc ) < 31). 



with little dependence on halo mass. If we assume that the gas frac- 
tions of these clusters are equal to the universal baryon fraction this 
implies that, for a 10^'' Mq cluster, the dust-to-gas ratio is ~ 
0.0003, or 3% of the interstellar medium ratio (Pei 199^. However, 
for a 10"'^ h'^ Mq group, the ratio is ~ 0.0055, approximately 
50% of the interstellar medium ratio. 

Here it is worth examining the theor etical expectations of ther- 
mal sputterine. lPraine & SalpeteJ i ll 9791) have shown that, for dust 
grains principally composed of graphite, silicate or iron, collisions 
with gas of 10^ < T < 10^ K can destroy them. Then, for gas of 
this temperature dust grains can have a typical lifetime r of 



2 X 10" yr 



0.01/im 



(6) 



where nn is the gas density and a is the radius of the typical dust 
grain. For the more massive clusters, a typical gas density is a 
few 10~'^ cm"*^ wh ile interstellar grains are smaller than 0.5 fim 
jMathis etal]|l977h . Therefore, dust grains in this gas will only 
survive for ~ 1 Gyr. So any dust which exists in the more massive 
clusters must have been accreted within the last Gyr. iMcGee et al.l 
l l2009h have shown, using semi-analytic merger trees, that local 
galaxy clusters accrete ~ 5% of their final galaxies per Gyr from 
otherwise isolated halos. This suggests that only 5% of the total ac- 
creted dust should exist in massive clusters today. This is remark- 
ably similar to the implied dust-to-gas ratio of the most massive 
clusters. 
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Figure 8. The excess number of galaxies (i < 21) as a function of cluster- 
centric position for three bins of halo mass. The excess number was deter- 
mined by statistical background subtraction. The uncertainties are 1 sigma 
expectations from the lackknife method. 



4.1 Excess galaxy profile 

iMenard et al.) l l2009bl) have recently shown that the dust-to-total 
mass ratio is approximately constant to a distance of 10 Mpc from 
individual galaxies. Similarly, they have shown that the dust-to- 
galaxy ratio is constant, except within the inner 20 kpc of the 
galaxy, where it's expected the signal is due to the galactic disk. 
We would like to explore similar ratios as a function of group cen- 
tric radii. Unfortunately, we do not have available weak lensing 
mass profiles for our groups and clusters. However, we are able to 
correlate the excess number of galaxies in the SDSS photometry 



catalogues with the positions of our groups and clusters. Addition- 
ally, it is worth pointing out that Sheldon et al. (2007a) has found 
that, for clusters similar to ours, the weak lensing mass to light ratio 
is approximately constant from 1 Mpc to ~ 20 Mpc. Therefore, we 
would expect that the galaxy/dust ratio should show similar trends 
to the mass/dust ratios. 

In Figure [8] we present the excess number of galaxies as a 
function of cluster-centric position for each of the three bins of halo 
mass. This measurement was made by correlating the position of 
all SDSS galaxies with a galactic de-reddened i magnitude of < 
21, where special care was taken with survey edges by using the 
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Figure 9. The excess surface dust mass per excess galaxy as a function of 
distance from the group centre. 



SDSSpix codeQ We choose to use the magnitude limit of i < 21 
because this limit was used by iMenard et al. (200_9b). This figure 
demonstrates that the excess number of galaxies has a steeper slope 
in the inner 10 Mpc than the similar excess dust mass plot. The 
figure also shows that there are significantly more galaxies in the 
higher halo mass bins, which again argues that the halo masses are 
well defined for large statistical samples. 

4.2 Dust mass associated with each galaxy 

Now that we have compiled the excess dust mass distribution 
around the clusters and the excess galaxy number around the clus- 
ters, we can examine their ratio. We show this in Figure|9] There are 
a number of striking aspects of this plot. In particular, we see that 
for the bulk of the galaxies, those which are >1 h^^ Mpc away, the 
excess dust mass per galaxy is on the order of lO' M©. This is ap- 
proximately equal to the stellar mass of the median excess galaxy, 
and roughly agrees with supernova feedback models, which imply 
that up to 50% of met gJs produced are blown out of the galaxy 
jFinlator & Davdl2008l) . The constant dust per galaxy fraction also 
confirms our earlier suggestion that the dust excess seen out to 
large radii (> 30 Mpc ) is simply due to excess halo clustering 
around galaxy groups and clusters. 

However, the most obvious feature of this plot is the signifi- 
cant drop in dust mass per galaxy at small groupcentric radii, espe- 
cially for the more massive cluster bin. The most massive clusters 
have dust to galaxy ratios ~ 15 % of the value far from the clus- 
ter, while small groups have ~ 80 % of the value. These numbers 
are not significantly different from the dust to gas ratios we found 
earlier. 

There are two explanations for this reduction in the dust to 
galaxy ratio. First, as we have discussed, the dust could be de- 
stroyed by thermal sputtering in hot gas. Within a fixed physical 
radius, this sputtering would be more effective in clusters than in 



http://dls.physics.ucdavis.edu/~scranton/SDSSPix/ 



groups, since the average density of the hot gas is substantially 
higher. Second, it is also possible that galaxies near clusters and 
groups may be fundamentally different in their dust creation prop- 
erties. It is expected that most of the dust creation occurs in su- 
pernova, and it has been shown through extensive observations 
that groups and clusters have significantl y lower star formation 
rates than isolated environments (Gomez et aI.ll2003l ; lBaIogh et al.l 
Pasquali et al. 2009) and that differ- 



2: 



2004: Wei nmann et alj|2006 ': 

ence extends to at l east z=l dWilman et al. 2005 : Ger ke et al. 20071 : 
iBalogh et alj2009l) . 

Unfortunately, the star formation rates in the large scale en- 
vironments around groups and clusters ha ve not been stud ied as 
extensively, especially at higher redshift. Lewis et al] ( 12002 ^ have 
attempted to quantify this in the low redshift universe using a sam- 
ple of massive clusters. They find that the mean star formation rate 
of galaxies is below the field value as far out as ~ 7-10 virial radii. 
Further, they see that the average star formation rate at 1 virial ra- 
dius is 70 % the isolated value. However, even despite this, the key 
to the amount of dust created is actually the stellar mass, which is 
essentially the integrated star formation rate. Since g alaxy groups 
are actually quite efficient at producing stellar mass jParker et al.l 
2005; Balogh et al. 2007), this suggests that the main driver of 
lower dust to galaxy ratios near clusters is not due to reduced dust 
creation. In effect, because the bulk of star formation/dust creation 
occurs at higher redshift the current cluster and field star formation 
rates are not as relevant. 



4.3 Impact on other science 

We have shown that there is a significant and large scale variation 
in the dust reddening associated with galaxy groups and clusters. 
Although studying the properties of this dust and its implications 
for dust creation and destruction properties was the main goal of 
the paper, in this section we briefly assess the impact this reddening 
and attenuation of l ight mi ght have on other studies. 

iMenard et alj ( l2009ah have shown that not correcting for the 
impact of dust on background "standard candle" supernova can bias 
the measurement of f^m at the few percent level. Here we point out 
an additional concern in the measuring of cosmological parameters. 
Our measurements imply that variations in attenuation by dust can 
be as high as 0.04 magnitudes in the V band (and higher in bluer 
bands). The next generation of large scale photometric surveys, 
such as LSST and DES, which principally are driven by the desire 
to probe the cosmological parameters, have design g oals of O.OI 
magnitude zeropoint stab ility from field to field (e.g. Ilvezic et"ail 
l2008l : lTuckeretal.ll2007l) . Our measurements suggest this will not 
be possible to achieve without accounting for the dust in and around 
groups and clusters. 

The dust attenuation may also be a problem for very low sur- 
face brightness features in clusters, suc h as the intra-cluster and 
-group medium. IZackrisson eTal] ( l2009h have shown that because 
the dust extinguishes background light, measuring the sky back- 
ground, of which the extragalactic background is a small part, 
some annulus from the cluster would lead to an over-subtraction. 
Crucially, this can lead to an underestimate of the intracluster 
light by as much as 1.3 magnitudes at an observed surface bright- 
ness of 30 magnitudes/arcsec'^ for attenuations of Av=0.Q5. 30 
magnitudes/arcsec^ is close to the faint lim it of the measurements 
of intracluster light by Izibetti et al.l ( |2005|) . While an attenuation 
of Av=Q-Q5 is not unreasonable given our results, as we have 
show n, the dust attenuat ion is dependent on the cluster-centric ra- 
dius. Izibetti et alj ^20051) et al. measured the sky background in a 
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lOO-kpc thick ring with an inner radius of 1 Mpc centered on their 
clusters. Although this is below the spatial resolution of our mea- 
surements, we do note that our most massive cluster varies only 
from Av ~ 0.042 at I h'^ Mpc to ~ 0.034 at 3 h'^ Mpc . In 
other words, we would expect that the variation in dust attenua- 
tions is not large (< 0.005 mag), and therefore does not induce a 
significant imderestimate in the intracluster Ught. 



5 CONCLUSIONS 

We have used a sample of spectroscopically identified high red- 
shift quasars to probe the dust content of 70,000 uniformly selected 
galaxy groups and clusters. We have used the resulting colour ex- 
cesses, in conjugation with synthetic dust models, to estimate the 
excess surface mass density aroimd these clusters. Finally we have 
compared the resulting dust distribution with the excess galaxies 
around the same groups and clusters. Our findings are as follows. 

• We have shown that there exists a large scale distribution of 
dust centered on groups and clusters with masses as low as M ~ 
10^^'^ Mq, and which extends 30 Mpc from the group 
centre. The wavelength-dependence of this extinction is consistent 
with that expected for interstellar dust. 

• We find that the dust must be distributed relatively uniformly, 
with a covering fraction on the order of 10 % to explain the excess 
color as a function of percentile of the color distribution. 

• We find that the halo mass dependence of the dust content is 
much smaller than would be expected by a simple scaling, imply- 
ing that the dust-to-gas ratio of the most massive clusters (~ 10^"* 

Mq) is ~ 3% of the local ISM value, while in small groups 
(~ lO^^-'' /i"^ Mq) it is ~ 55% of the local ISM value. 

• We find that the implied dust-to-galaxy ratio falls significantly 
closer to the group and cluster center. This reduction in the dust to 
galaxy ratio has a significant halo mass dependence, such that the 
more massive groups and clusters show a stronger reduction. This 
suggests that either dust is destroyed by thermal sputtering of the 
dust grains by the hot, dense gas or the intrinsic dust production is 
reduced in these galaxies. 
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